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Abstract 

We study the long-time behavior of decoupled continuous-time random walks characterized by 
superheavy-tailed distributions of waiting times and symmetric heavy-tailed distributions of jump 
lengths. Our main quantity of interest is the limiting probability density of the position of the 
walker multiplied by a scaling function of time. We show that the probability density of the scaled 
walker position converges in the long-time limit to a non-degenerate one only if the scaling function 
behaves in a certain way. This function as well as the limiting probability density are determined 
in explicit form. Also, we express the limiting probability density which has heavy tails in terms 
of the Fox i?-function and find its behavior for small and large distances. 
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I. INTRODUCTION 



Continuous-time random walks (CTRWs), introduced by MontroU and Weiss [1], con- 
stitute an important class of jump processes that are widely used to model a variety of 
physical, geological, biological, economic and other phenomena. In particular, these pro- 
cesses describe anomalous diffusion and transport in disordered media (see, e.g., Refs. [2-4] 
and references therein), seismic [5, 6] and financial [7, 8] data. A remarkable fact is that sys- 
tems so different from one another can successfully be described within the CTRW approach. 
This is because two random variables that many systems have in common, the waiting time 
between successive jumps and the jump length, are used to model the CTRW. Therefore, 
even the decoupled CTRW, when these variables are independent, is rather flexible. 

The probability density P{x, t) of the walker position X{t) is the most important charac- 
teristic of the CTRW. It satisfies the integral master equation [9-11] which in the decoupled 
case depends only on the probability density p(r) of waiting times and on the probability 
density w{x) of jump lengths. Because exact solutions of this equation are known in very 
few cases [12-15], there is considerable interest in studying the long-time behavior of P{x, t) 
that is responsible for the transport and diffusion properties of objects described by the 
CTRW model. In this context, much attention has been paid to the probabihty densities 
p{r) and w{x) having finite second moments and/or to those having heavy tails. It has been 
established [16-18] that different combinations of these properties of the waiting time and 
jump densities lead to different long-time distributions of X{t). In Ref. [19], all possible 
distributions were expressed in terms of the limiting distributions of the properly scaled 
walker position. 

In some cases the waiting-time densities are assumed to be superheavy-tailed, i.e., such 
that all fractional moments of p(r) are infinite. In particular, this class of densities is used 
to model the superslow diffusion in which the diffusion front spreads more slowly than any 
positive power of time [20-23]. In general, one might expect superheavy-tailed distributions 
to reflect extremely slow time-dependent phenomena such as may occur in some relaxation 
and aging processes. Such distributions are also applicable within a Langevin rather than a 
CTRW description when dealing with processes that are interrupted by an absorption event 
or by the transition of a particle to a qualitatively different state [24, 25]. The long-time 
behavior of the decoupled CTRWs characterized by these waiting-time densities and jump 
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densities with finite second moments is considered in Ref . [26] . Here we focus on asymptotic 
solutions of the CTRWs in the case when the densities p(t) and w{x) are superheavy- and 
heavy-tailed, respectively. 

The paper is structured as follows. In Sec. 11, we formulate the main definitions and write 
the basic equations describing the decoupled CTRW. A one-parameter limiting probability 
density of the scaled walker position that corresponds to the superheavy-tailed distributions 
of waiting times and the symmetric heavy-tailed distributions of jump lengths is determined 
in Sec. III. Here, we also find the scaling function and prove the positivity and unimodality 
of the limiting probability density. In Sec. IV, we express the limiting density and the 
corresponding cumulative distribution function in terms of Fox iJ-functions and consider a 
few particular examples. The short- and long-distance behavior of the limiting density is 
studied in Sec. V. Our main results are summarized in Sec. VI. 



II. MAIN DEFINITIONS AND BASIC EQUATIONS 

The CTRW approach deals with a wide class of continuous-time jump processes X{t) 
represented as 

Nit) 

= (2.1) 

n=l 

Here, N{t) = 0, 1, 2, ... is the random number of jumps that a walker has performed up to 
the time t (if N{t) = then X{t) = 0), and Xn G (— cxo, cxd) are the independent random 
variables (jump lengths) distributed with some probability density w{x). In order to specify 
the counting process N{t), the waiting times t„, i.e., times between successive jumps, are 
introduced. Like the jump lengths, the waiting times are assumed to be independent random 
variables distributed with probability density p(t). If the variables x„ and t„ are independent 
of each other as well, i.e., if the CTRW is decoupled, then the probability density P{x, t) of 
the walker position X{t) depends only on w{x) andp(r). According to [1], in Fourier-Laplace 
space this dependence has the form 

Pks = y (2.2) 

s{l-psWk) 

where Wk — J-{w{x)} — J^^dxe'^'^^w{x) (— oo < k < oo) is the Fourier transform of 
w{x), Ps = C{p{t)} = dte-'*p{t) (Res > 0) is the Laplace transform of p(t), and 
Pks^J'{C{P{x,t)}}. 
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Prom Eq. (2.2) one can get 

Ps{x) = ^^-^ I I + Six) (2.3) 

S [1- PsWk J S 

and 

P{x, t) = f (l-P^)P^ ^-i| ^^f^ II + ^(^)^(^) _ (2.4) 

Is I 1 - Psti^fe J J 

Here, J^~^{fk} = fix) = (27r)^^ (i/ce~*^^/fe is the inverse Fourier transform, 6{x) is the 
Dirac S function, C~^{gs} = g(t) = (27ri)~^ J^^^^ dse^^Qg (c is a real number exceeding the 
real parts of all singularities of gs) is the inverse Laplace transform, and 

V(t) = C-'l^^^^ = J^^ drpir) (2.5) 

with V{0) = 1 and V{oo) = is the survival or exceedance probability. Using Eq. (2.4), 
the integral formula dxe^^^^ = 2n6{k) and the well-known properties of the 5 func- 
tion, it is not difficult to show that the probability density P{x, t) is properly normalized: 
J^^dxP{x,t) — 1. Since X{0) — 0, the initial condition for P{x,t) reads P{x,Q) — 5{x) 
and, if boundary conditions are not imposed, P{x, i) — > as i — > oo. 

According to this last property, the probability density of the walker position vanishes 
in the long-time limit. It is therefore reasonable to introduce the scaled walker position 
Y{t) = a{t)X{t) and find the positive scaling function a{t) such that the limiting probability 
density 

P(y) = lim 4t^( (2-6) 

of Y{t), i.e., the probability density of the random variable Y[oo), is non- vanishing and 
non-degenerate. The importance of the functions a{t) and V{y) is that, since P{x,t) ~ 
a{t)V{a{t)x) as t — )■ oo, they completely describe the long-time behavior of the original 
walker position X{t). To satisfy the above requirements on V{y), the scaling function must 
go to zero as f — > oo in a certain way. In fact, these requirements permit one to determine a{t) 
up to a constant factor which, however, is not important and can be chosen for convenience. 

The pairs a{t) and V{y) have been determined for all cases characterized by finite second 
moments and/or heavy tails of the probability densities p(r) and w{x) [19]. In contrast, 
the case with superheavy tails has been much less studied. In fact, the pair a(t) and V{y) 
has been determined only when p{t) has a superheavy tail and w{x) has a finite second 
moment I2 [26]. Because I2 — 00 ii w{x) is heavy tailed, one may expect that in this case 
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the long-time behavior of the walker position changes qualitatively and thus the pair a{t) 
and 'P(y) changes as well. More precisely, in this paper we study the long-time behavior of 
decoupled CTRWs whose waiting-time densities p(r) and jump densities w{x) [it is assumed 
that w{—x) = w{x)] are described by the asymptotic formulas 

(r ^ oo) (2.7) 

(kl^oo), (2.8) 

where the positive function h{T) varies slowly at infinity, i.e., h{^T) ~ h{T) as r — >■ oo for 
all // > 0, the tail index a is restricted to the interval (0, 2], and w > 0. The waiting-time 
and jump densities considered here belong to the classes of superheavy- and heavy-tailed 
densities, respectively. The difference between these classes consists in different asymptotic 
behavior of the constituent probability densities that, in turn, results in different properties 
of their fractional moments. Specifically, while the fractional moments drr^p^r) of p(r) 
are infinite for all p > 0, the fractional moments J^^dx\x\'^w{x) of w{x) are infinite only 
ii p > a. It should also be noted that the conditions u > and a e (0, 2] are completely 
compatible with the normalization condition J^^dxw{x) — 1. In contrast, the normalization 
condition dTp{T) = 1 imposes an additional restriction on the asymptotic behavior of 
/i(r): /i(r) = o(l/ Inr) as r — )■ oo. 



PKV 

T 

and 

u 



III. SCALING FUNCTIONS AND THE LIMITING PROBABILITY DENSITY 

According to the Tauberian theorem for Laplace transforms [27], the long-time behavior 
of the probability density P{x,t) is determined by the asymptotic behavior of the Laplace 
transform Ps{x) when the real parameter s tends to zero. Because the waiting-time distri- 
bution is normalized to unity, the condition — > 1 holds as s — >■ 0. It follows from Eq. (2.3) 
that wc also need to find the s — )■ behavior of 1 — Ps- To this end, it is convenient to 
use the representation 1 — ps = Jl^ dqe~W {q / s) which, together with the fact [26] that the 
survival probability V{t) varies slowly at infinity, immediately gives 

l-Ps^V{l/s) (3.1) 
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as s — > 0. Then, taking into account that as s — > the main contribution to T — 
Vs^k)) comes from a small vicinity of the point A; = 0, i.e., 

\\-l.w}\ ~ \v{\ls) + l-w\ ^^-^^ 
Eq. (2.3) in the small-s limit yields 

The long-time behavior of P{x^ t) can be found directly from the limiting formula (3.3) 
by applying the above mentioned Tauberian theorem. It states that if the function v{t) is 
ultimately monotonic and Vg ~ s~'^L{l/s) (0 < 7 < 00) as s — >■ 0, then v{t) ~ t'^~^L{t) /r{'y) 
as t — >■ 00. Here, r(7) denotes the gamma function and L{t) is a slowly varying function at 
infinity. In our case 7 = 1, therefore from Eq. (3.3) one obtains 

X xc-/ X Vii) f°° „ cos(xA;) 

Jo V{t) + l-Wk 

{t 00). Since in the long-time limit (when V{t) tends to zero) the main contribution to 
the integral in Eq. (3.4) comes from a small vicinity of the point k — 0, the exact formula 

POO 

1 — Wk = 2 / dx[l — cos(A;a;)]'u;(a;) (3.5) 
Jo 

can be replaced by one valid in this regime. Using Eq. (2.8) and the integral relation 

1-C0S{X) _ TT 

Jo x^+'^ 2r(l + a) sin(7ra/2) ^ ' ^ 

{0 < a < 2), from Eq. (3.5) at \k\ ^ we find 

TT?/ 

r(l -I- a) sm(7ra/2j 

Substituting this result into the asymptotic formula (3.4) and applying the definition 
(2.6), the limiting probability density V{y) can be written in the form 

Viy) = Um i r d. ^ -. (3.8) 

^ r(l+a)sm(7ra/2)y(t) 

It appears from this that V{y) is non-vanishing and non-degenerate only if the factor in 
front of tends to a nonzero finite limit as i — > 00. Assuming for convenience that this 
limit equals 1, we obtain the asymptotic representation of the scaling function 



{t — > oo) and the corresponding limiting density 

Viy) = ^rdx'-^ (3.10) 

(the fact that V{y) is a probabihty density will be proved below). The symmetry condition 
'Pi~y) — T^{y)-i which follows from Eq. (3.10), is a consequence of the symmetry of the jump 
density w{x). 

Since at a = 2 the integral in Eq. (3.6) diverges, the limiting formula (3.7) is not applicable 
to this case. Therefore, in order to find 1 — at a = 2 and \k\ — > 0, wc first split the 
interval of integration in Eq. (3.5) into two parts, (0, h) and (6, oo) with 6 ~ 1. Then, taking 
into account that as — >■ the contribution of the first interval to the right-hand side of 
Eq. (3.5) can be approximated by k"^ dxx'^w(x) and the second one by uk'^ln.(l/\k\), we 
get 

1 - ujfc ~ In ^ (3.11) 
\k\ 

(|A;| ^ 0). In accordance with this, the limiting probability density when a — 2 takes the 
form 

V(y) - lim - r dx . (3.12) 

As before, we choose the long-time limit of the factor in front of to be equal to unity. 
In this case the asymptotic behavior of the scaling function a(i) is determined by the relation 



ua}{£) ln[l/a(i)] ~ y(t) {t — >■ oo). Assuming that a{i) ~ ^JV{t)|ua\{t\ where the positive 



function ai(t) satisfies the conditions a\{i) — >■ and yJVit) — o(ai(t)) as i — >■ oo, from this 
relation we obtain a^it) ~ -^2/ ln[l/1/(t)], and thus 



{t oo). The limiting probability density (3.12) which corresponds to this scale function 
is given by 



V{y) = 3 / dx^-^^ = ^e-l^ (3.14) 



TT io l + X^ 2 

showing that Eq. (3.10) is valid for a = 2 as well. Wc note that the same two-sided 
exponential density (3.14) describes the limiting distribution when the jump density w{x) 
has a finite second moment I2 — J^^dxx^w{x) [26]. However, because at Z2 < 00 the 
asymptotic behavior of the scaling function, a(t) ~ ■\/2V{t) /I2, is quite different from that 
given in Eq. (3.13), the long-time behaviors of the walker position in these cases are also 
quite different. 
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A. Positivity and unimodality of 'P{y) 



To be a probability density, the function 7'(y) must be normalized and positive (non- 
negative). The normalization condition J^^dyV{y) — 1 can easily be proved using 
Eq. (3.10), which represents V{y) as a cosine Fourier transform, and the integral repre- 
sentation S{x) = (l/27r) J^^dycos{yx) of the S function. However, except for the case 
a = 2, where according to Eq. (3.14) V{y) > 0, the use of Eq. (3.10) to prove the positivity 
of Vly) is impractical because of the oscillating character of the integrand. On this point, 
the representation of V{y) in the form of a Laplace transform would be preferable. In order 
to find it, we first define the function 

1 pi\y\z 

m - -V—^ (3-15) 

(0 < q; < 2) of the complex variable z — x+iu. This function is analytic in the first quadrant 
of the z-plane (when l-zj > and < arg2; < 7r/2), and so from the Cauchy integral theorem 
[28] we have dzf{z) — 0, where C is a simple closed contour that lies in the domain of 
analyticity of f{z). Then, choosing the contour C to be the boundary of the first quadrant 
(we emphasize that the branch point 2; = is outside the contour) and applying the Jordan 
lemma [28], the above integral reduces to 

poo roG 

/ dxf{x)-i duf{iu)^0. (3.16) 
Jo Jo 

Finally, taking into account that V{y) — Re [J^^ dxf{x)] and = cos(7rQ;/2) -|- i sin(7ra/2), 
from the real part of Eq. (3.16) we obtain 

TT Jo l + 2cos(7ra/2)x" + x2" ^ ^ 

The main advantage of this representation of V{y) is that it clearly shows that V{y) > 
when < q; < 2. Thus, since V{y) is positive for a = 2 as well, we can conclude that 
the function V{y) is indeed the probability density for all a in the interval (0,2]. Another 
important property of V{y)^ which follows directly from Eq. (3.17), is that dV{y)/dy < 
when y > 0. Together with the condition V{—y) = V{y), it shows that the limiting 
probability density is symmetric, unimodal and centered at the origin. In contrast to the 
scaling function, which depends on all the parameters characterizing the asymptotic behavior 
of the waiting time and jump densities, the limiting density depends only on the tail index 
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FIG. 1: Plots of the probability density V{y) for two values of the tail index a belonging to the 
intervals (0,1] and (1,2]. 

a. According to Eqs. (3.10) and (3.17), this parameter strongly influences the properties of 
V{y). In particular, the behaviors of V{y) in the vicinity of the origin differ substantially 
from one another when a G (0,1] and a G (1,2], as illustrated in Fig. 1 (for details, see 
Sec. V). 



IV. LIMITING DISTRIBUTION IN TERMS OF SPECIAL FUNCTIONS 

To get more insight into the mathematical structure of the limiting probability density 
V{y), it is reasonable to express it in terms of well-known special functions. Toward this 
end, we flrst represent 'P(y) as the inverse Mellin transform. The Mellin transform of a 
function f{y) is deflned by fr = M.{f{y)} = dyf{y)y'^^^. Therefore, for the function 
f{y) = dxu{yx)v{x) one gets fr = Uj-Vi^r [29]. If f{y) is associated with V{y) from 
Eq. (3.10), then the functions u{x) and v{x) can be chosen as u{x) = vr"^ cos(x) and v{x) = 
(1 -|- x")~^ whose Mellin transforms are given by [29] 

Ur = -r(r) cos(^yj (0 < Rer < 1) (4.1) 

and 



= -y(L]y(i--] (0<Rer<a). (4.2) 

Using the reflection formula [30] r(l/2— r/2)r(l/2-|-r/2) = vr/ cos(7rr/2) to replace cos(7rr/2) 

in Eq. (4.1), the Mellin transform Vr = Ur-Vi^r of 'P{y) takes the form 

^ r(r)r(l-l/a + r/a)r(l/a-r/a) 

ar(l/2-r/2)r(l/2 + r/2) ' ^ ' 



where max(l — q;,0) < Rer < 1. Finally, introducing the inverse Mellin transform as 
M-'{fr} = f(y) = (27ri)-i dr/.y-'' and utilizing the fact that V{-y) = V(y), we find 



rc+ioo 



(4.4) 



The structure of Vr suggests that the probabihty density V{y) is a particular case of the 
Fox if-function which can be defined by means of a Mellin-Barnes integral as follows (see, 
e.g., Ref. [31]): 



Here, 



Tjm.n 



{ai,A 



1), - ■ ■ , 



(6i,5i),...,(6„i?,) 

n7=ir(6, + fi,r)n;ur(i 



(4.5) 



(4.6) 



h, - B.r) , . r(a,- + Af) ' 
m, g are whole numbers, < m < g, < n < p, and 6j are real or complex numbers, 
Aj, Bj > 0, and L is a suitable contour in the complex r-plane which separates the poles of 
the gamma functions r{bj + Bjr) from the poles of the gamma functions r(l — aj — Ajr). 
It is also assumed that an empty product equals 1. Comparing Eqs. (4.3) and (4.4) with 
Eqs. (4.6) and (4.5), respectively, we see that 



a 



{l-l/a,l/a), (1/2,1/2) 
(0,l),(l-l/a,l/a),(l/2,l/2) 



(4.7) 



It should be noted that the cumulative distribution function F{y) — 1/2 + dy'V{y') 
of the random variable y(oo) can also be expressed through the if-function. To show this, 
we write F{y) — dy'V{y') and take into account the following property of the Mellin 
transform [29]: M {Jq f{y')dy'} = —fr+i/r. According to this, Fr = —Vr+i/r and, from 
Eq. (4.3) and the functional equation r(l + ,x) = xr{x), one gets 

r(r)r(r/a)r(l -r/a) 



ar(l -r/2)r(r/2) 

[— min(l,Q;) < Rer < 0]. Therefore, using Eqs. (4.5) and (4.6), we obtain 



(4.8) 



1 
2 



sgn(y) 



H. 



a 



2,1 
2,3 



(0,1/a), (0,1/2) 
(0,1),(0,1/q;), (0,1/2) 



(4.9) 



10 



A. Particular examples 



For some special values of the tail parameter a the if-functions in Eqs. (4.7) and (4.9) can 
be reduced to more familiar special (or even elementary) functions. Because the probability 
density V{y) and the distribution function F{y) provide equivalent descriptions of the long- 
time behavior of the scaled walker position Y{t), next we consider only the properties of 
Viy). The simplest situation occurs when a = 2. In this case both reduction formulas [31] 
can be applied, yielding 



1 ^^2,1 



1 0-2,0 



1,0 



(1/2, 1/2), (1/2, 1/2) 
(0,1), (1/2, 1/2), (1/2, 1/2) 

(1/2,1/2) 
(0,1), (1/2, 1/2) 



(0,1) 



(4.10) 



Since the last if-function equals e 1^' [31], this ascertains that Eq. (4.7) at a = 2 reduces to 
Eq. (3.14). 

If the parameter a is rational, then the probability density Viy) can, in principle, be 

expressed in terms of the Mcijcr G-function as well. The G- function, which is a particular 
case of the i7-function, is defined as 



p,q 



Oi, . . . , Clp 
bi,...,bg 



27rz Jl 



(4.11) 



with — ^r\Aj,Bj=i- As a first illustrative example, we consider the case when a — 
1. Changing the variable of integration in Eq. (4.4) from r to 2r, one readily obtains 
V{y) = (ni)-^ drV2r{y^)~% where max(l - a,0)/2 < c < 1/2. Then, using Eq. (4.3) 

with a = 1 and the duplication formula [30] T{2r) = 7i-^/'^2^'-^T {r)T {1/2 + r), the Mellin 
transform V2r can be written in the form 

22r 



2r 



47^3/2 



r^(r)r(l/2 + r)r(l -r). 



(4.12) 



Therefore, in accordance with the definition (4.11), the limiting probability density (4.7) at 
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a — 1 has the following G-function representation: 



2,1 



(0,1), (1/2, 1/2) 
(0,1), (0,1), (1/2, 1/2) 



27r3/2 



^1,3 



I/V4 





0,0,1/2 



(4.13) 



Remarkably, the limiting probability density V{y) a.t a = 1 can be expressed not only in 

terms of the Fox and Meijer functions, but also in terms of the well-known sine, si(y) = 
— dx sm{x) / X , and cosine, Ci{y) = — dx cos{x)/x, integral functions. Indeed, using 
the exact result for the cosine Fourier transform of (1 + x)~^ [32], we obtain 



r{y) = — [sin(|y|)si(||/|) + cos(y) Ci{\y\)]. 



TT 



(4.14) 



Finally, in our last example we consider the case a = 1/2. Following straightforward 
calculations similar to those described above, for the Mellin transform V2r we obtain the 
expression 

22r 



-r(-l/4 + r)r'(r)r(l/4 + r)r(l/2 + r) 



xr(5/4 - r)r(l - r)r(3/4 - r) 



(4.15) 



from which it follows that 



V{y)^2H, 



2,1 
2,3 



(-1,2), (1/2, 1/2) 
(0,1), (-1,2), (1/2, 1/2) 



47r7/2 



G 



5,3 
3,5 



yV4 



-1/4,0,1/4 
-1/4,0,0,1/4,1/2 



(4.16) 



V. ASYMPTOTIC BEHAVIOR OF V{y) 

Using Eq. (4.7), the behavior of the limiting probability density 'P{y) for small and large 
values of |y| can, in principle, be found from the expansions obtained for the i7- function 
in different hmits (for details, see Ref. [31] and references therein). However, because V{y) 
is a very particular case of the if-f unction, it is reasonable and convenient to derive the 
corresponding limiting formulas directly from the source representation (3.10). 
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A. Short-distance behavior 



There are three regions of the tail index a, which we consider separately, where the 
limiting behaviors oiV{y) as \y\ — > differ from one another. 

q; e (0, i). In this case Eq. (3.10), after changing the variable of integration from x to 
as \y\ yields 

^/ N 1 /"°° , cos(x) 

Then, since Jq°° dxcos{x)/x°' — T{1 — a) sin(7rQ;/2), one gets 

^ r(l-.)sin(W2) ^ ,5 2) 

'^\y\ 

a = 1. Using the formulas si(|?/|) ~ \y\ and Ci(|y|) ~ ln|7/| (|y| ^ 0) [33], Eq. (4.14) 
which follows from Eq. (3.10) immediately yields 

V{y) ^ --\Ti\y\. (5.3) 

TT 

a & {1,2\. Finally, in this case it is convenient to rewrite Eq. (3.10) in the form 

W = nO)-^ rdx'-f^^, (5.4) 

where 'P(O) = [a sin(7r/a)]^^. Then, neglecting \y\" in the integrand and taking into account 
that dx[l — cos{x)]/x°' — r(2 — a) sm{iTa/2)/{a — 1), we obtain 

Pfa)^_L^- r(2-«)sin(W2) 

Q;sm(7r/Q;) 7r(Q; — 1) 

It should be noted that, since lima;_).o r(x) sin(7rx/2) = 7r/2, the limiting formula (5.5) at 
a = 2 reduces to V{y) ~ (1 — \y\)/2, in accordance with Eq. (3.14). 



B. Long-distance behavior 

The asymptotic behavior of V{y) as |y| — >■ oo can easily be found by a single (if < a < 1) 
or double (if 1 < a < 2) integration by parts of Eq. (3.10) with a subsequent change of the 
integration variable from x to x/\y\. In particular, for a e (0, 1) this yields 



7r||/|i+°7o x^-'-[l + {x/\y\Y 



V fdx^ (5.6) 



T^\y\ Jo ^ 
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and so 

^ r(l + a) sin(7ra/2) 

It is not difficult to verify that the asymptotic formula (5.7) also holds for a G (1,2). 
Moreover, since Eq. (4.14) leads to V{y) ~ 7r~^|y|~^ as \y\ oo, this formula is valid for 
q; = 1 as well. 

Thus, according to Eq. (5.7), the hmiting probabihty density V{y) when a e (0,2) is 
heavy-tailed with the same tail index a as in the jump density w^x). In contrast, at a = 2 
the limiting density has exponential tails, while the jump density is still heavy-tailed, see 
Eq. (2.8). We also note that the same tail index a characterizes the limiting probability 
density when both the waiting-time and jump distributions are heavy-tailed [19]. However, 
this does not mean that the long-time behaviors of the CTRWs with heavy- and superheavy- 
tailed distributions of waiting times are identical. This is because the scaling functions for 
these CTRWs are quite different. Specifically, while in the former case the scaling functions 
are power functions of time [19], in the latter case they vary more slowly, see Eqs. (3.9) and 
(3.13). 



VI. CONCLUSIONS 



We have determined a new class of asymptotic solutions of the CTRWs characterized by 
superheavy-tailed distributions of waiting times and symmetric heavy-tailed distributions of 
jump lengths. These solutions represent the probability densities of the scaled walker posi- 
tion, i.e., the random walker position multiplied by a time-dependent deterministic scaling 
function, in the long-time limit. We have found both the limiting probability densities and 
the corresponding scaling functions which completely describe the long-time behavior of the 
reference CTRWs. It turns out that the scahng functions depend on the survival probabil- 
ity characterizing the long-time behavior of the waiting-time density and on the tail index 
a e (0, 2] describing the asymptotic behavior of the jump density. In contrast, the limit- 
ing densities, which have been represented in the form of Fourier and Laplace transforms, 
depend only on a. 

The limiting probability densities Viy) form a class of symmetric and unimodal functions 
centered at the origin. Among other things, we have determined the limiting behavior of 
these densities for small and large distances. We find that while at a = 2 the function V{y) 
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has exponential tails, at o; e (0, 2) the tails are heavy and are characterized by the same tail 
index a as the jump density. In the vicinity of the origin, the behavior of V{y) for a e (0, 1] 
is quite different from that for a e (1, 2]. Specifically, P(0) is infinite in the former case and 
is finite in the latter. Finally, we have expressed the limiting probability densities in terms 
of the Fox if- function for the general case of arbitrary a and, for a few values of a, in terms 
of the Meijer G-function. 
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